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Abstract 

We present a probabilistic approach for the study of systems with exclusions, in the regime 
traditionally studied via cluster-expansion methods. In this paper we focus on its application 
for the gases of Peierls contours found in the study of the Ising model at low temperatures, but 
most of the results are general. We realize the equilibrium measure as the invariant measure of 
a loss-network process whose existence is ensured by a subcriticality condition of a dominant 
branching process. In this regime, the approach yields, besides existence and uniqueness of 
the measure, properties such as exponential space convergence and mixing, and a central limit 
theorem. The loss network converges exponentially fast to the equilibrium measure, without 
metastable traps. This convergence is faster at low temperatures, where it leads to the proof 
of an asymptotic Poisson distribution of contours. Our results on the mixing properties of the 
measure are comparable to those obtained with "duplicated- variables expansion" , used to treat 
systems with disorder and coupled map lattices. It works in a larger region of validity than 
usual cluster-expansion formalisms, and it is not tied to the analyticity of the pressure. In fact, 
it does not lead to any kind of expansion for the latter, and the properties of the equilibrium 
measure are obtained without resorting to combinatorial or complex analysis techniques. 
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1 Introduction 



In this paper we develop a probabilistic approach to the study of the equilibrium measure of 
systems with exclusions — like hard-core gases, contours, polymers or animals — in the low- 
density or extreme-temperature regime. This regime has traditionally been studied via cluster- 
expansion methods, which relied either on sophisticated combinatorial estimations [Malyshev 
(1980), Seller (1982), Brydges (1984)] or on astute inductive hypotheses plus complex analysis 
[Kotecky and Preiss (1986), Dobrushin (1996, 1996a)]. 

In contrast, we realize the equilibrium measure as the invariant measure of a loss network 
process that can be studied using standards tools and notions from probabilistic models and 
processes. Loss networks, first introduced by Erlang in 1917, encompass a rather general family 
of processes as discussed in Kelly (1991) and references therein. Technically we work with the 
so-called fixed-routing loss networks. We build on ideas of Ferrari and Garcia (1998) to show 
(Section ^^) that the existence of the loss network can be related to the absence of percolation 
in an oriented percolation process. This condition also yields other properties of the process 
and its invariant measure, like uniqueness, convergence through sequences of finite volumes 
and mixing properties (Theorem |4.6| ). More precise results can be obtained by resorting to 
a dominant multitype branching process (Section |5|). Roughly speaking, the mean number of 
branches of this process becomes the driving parameter: Subcriticality is a sufficient condition 
for the construction to work. Time and space rates of convergence and mixing rates are explicitly 
obtained in terms of this parameter (Theorem p. 101 ). 

The approach of this paper was already exploited in Fernandez, Ferrari and Garcia (1998). 
Here we refine and complete the theory presented there, extending their region of validity and 
including proofs of exponential mixing and convergence to a Poisson process. 

For concreteness, we analyze in this paper the gases of Peierls contours used, for instance, for 
low-temperature studies of the Ising model [Peierls (1936), Dobrushin (1965), Griffiths (1964)]. 
The subcriticality condition of the corresponding branching process is 
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where I7I indicates the length (perimeter, surface area) of a contour 7, w{'y) is its weight and 
"~" stands for the volume-exclusion (=non-intersection) condition. 

Condition (|1.1D is considerably weaker than those obtained by the most developed cluster- 
expansion approaches [Kotecky and Preiss (1986), Dobrushin (1996, 1996a)] in which the factor 
1^1 in the right-hand side is replaced by a function that grows exponentially with |^|. The weak- 
ening of the condition has a price: unlike previous approaches, ours does not yield analyticity 



properties of the expectations. Condition (1.1) is similar to conditions obtained in the study of 
systems for which nonanalyticity is known [von Dreifus, Klein and Perez (1995)] or suspected 
[Bricmont and Kupiainen (1996, 1997)]. Our results on space convergence and mixing rates 
can basically be obtained with the "duplicate-system" expansions of Bricmont and Kupiainen 
(1996, 1997). 

The novelty of our approach lies in the following features: First, it offers a completely differ- 
ent framework for the study of hard-core measures, based on well known stochastic processes. 
This can conceivably lead to new insights and stronger results. In particular the condition ( |l.lD 
can potentially be weakened via sub-criticality estimates obtained directly for the associated 
oriented percolation process, without resorting to a dominating branching process. Similar im- 
provements have been done for the contact process and oriented percolation, see Liggett (1995), 
for instance. Second, our construction involves a stochastic process (the loss network) that 
converges exponentially fast to the sought-after measure. This process is, in principle, easy to 
simulate and its potential as a computational tool deserves to be explored (Fernandez, Ferrari 
and Garcia, 1999). Its rate of convergence to the equilibrium measure increases as the tempera- 
ture decreases and, unlike spin-flip dynamics, it does not present meta-stable traps (any contour 
lives an exponential time of mean one). Finally, the construction permits a rather straightfor- 
ward proof of the asymptotic Poisson distribution of contours at low temperature (Theorem 



2.17 ). It is not obvious to us how such a result can be obtained through the standard statistical 



mechanical expansions. 

The present approach does not lead to a series expansion for the pressure (or free-energy 
density). In particular, it does not yield the "surface-tension bounds" that play such a crucial 



3 



role in some applications of cluster expansions [see eg. Zahradm'ck (1984), Borgs and Imbrie 
(1989)]. In fact, the approach is designed so to bypass expansions of this type. It is a probabilistic 
approach designed to answer probabilistic questions — existence of expectations, properties of 
correlation functions — in a direct way, without combinatorial or complex-analysis techniques. 
It is not an alternative to cluster expansions: it has a different regime of validity and different 
aims. 

While the ideas behind our results are natural and simple, their formalization requires many 
intermediate technical results that may obscure the development of our theory. Let us, therefore, 
present a sort of "road map" of the paper to guide the reader. The main results are presented, in 
a self-contained manner, in Section ^. The actual construction of the loss network is the subject 
of Section |3|. We start with a reference "free" process of Poissonian births and exponentially 
distributed deaths, with respect to which the loss network is absolutely continuous. The novelty 
of our approach resides in the fact that, rather than independently generating birth and death- 
times, the lifetimes are associated to each birth-time as a mark. Hence, unlike comparable 
constructions, each death-time has an associated birth-time and defines an space-time "cylinder" 
representing the presence of a loss event (contour). This permits to comb the process either 
backwards or forward in time with equal ease. The "time backwards" point of view leads to 
the notion of "backwards oriented percolation" which will be our main conceptual and practical 
tool. The idea is to construct the loss network by erasing from the free process those cylinders 
that conflict with pre-existing ones. This can only be done if the set of preexisting cylinders (the 
"clan of ancestors" ) is finite with probability one. This is precisely the condition of absence of 
backwards oriented percolation. The main technical result of this Section |3| is the proof that this 



construction yields precisely the loss network. [This is contained in the proof of Theorem 3.20|. 
It is apparent that the same construction works for any other process absolutely continuous with 
respect to a Poisson birth and death process, for instance point or Boolean processes (Baddeley 
and van Lieshout, 1995; Kendall, 1997, 1998). 

Once the process is so constructed, the course of action is clear. First, we relate time and 
space mixing properties with the time and space size of the percolation clan. This is done in 
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Section [Theorem 4^ ] . The only shghtly involved part of this section is related to the proof 
of the mixing properties. Indeed, we resort to a — standard but unusual — coupling between 
clans, together with a continuous-time construction (Section |4.5D , to improve a little over results 
proved previously by the method of "duplicated variables" (von Dreifus, Klein and Perez, 1995; 
Bricmont and Kupiainen, 1996). 

In order to present quantitative estimates in terms of the parameters of the problem, we fol- 
low the known technique (used, for instance by Hall, 1985) of bounding percolation probabilities 
via a branching process. In particular, the size of the percolation clan is bounded by number 
of branches. This is done in Section ^, where the only (very simple!) algebraic calculations of 
the paper are presented [displays ( p. 15] ), ( 5.41 ) and ( 5.44D ]. As expected, subcriticality of the 
branching process implies lack of percolation and exponentially damped sizes of the percolation 
clan. The main estimations are contained in Theorem ^.27 . Once again, we present a slightly 
unusual continuous-time construction (Section |5.2| ) to improve one of the estimations, namely 
the time- length of a clan [part (ii) of Theorem ^.27| . Readers can opt instead for the more 
direct, but slightly weaker, estimate presented in the first remark following the theorem. 

The results of Theorem 2.10 are a direct consequence of the estimates of Section |5| applied 
to the percolation expressions of Section ^, as explained in Sections |^ and |^ The proof of the 



Poisson approximation of the loss network (Theorem 2.17) requires some further considerations 
presented in the final Section |8|. 

Our work was motivated in part by the posthumous review of Roland Dobrushin (1996a), 
where he complained that while "perturbation methods are intensively used by mathematical 
physicists, they are not so popular as correlation inequalities among the probabilists" . He called 
for a "systematical exposition oriented to the mathematicians" . In this paper we follow his call: 
we concentrate on probabilistic issues and exploit probabilistic arguments. 



2 Definitions and results 
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2.1 The contour model 

We consider the d-dimensional lattice Z*^ and call plaquettes the {d — l)-dimensional unit cubes 
centered at points of Z'^. We identify each plaquette with its center. A set of plaquettes is called 
surface. Two plaquettes are adjacent if they share a (d — 2)-dimensional face. This defines a 
notion of connection: a set 7 of plaquettes forms a connected (hyper) surface if for every two 
plaquettes x,y there is a sequence of pairwise adjacent plaquettes starting at x and ending at y. 
A surface is closed if every {d — 2)-dimensional face is shared by an even number of plaquettes 
in the surface. A contour 7 is a connected and closed family of plaquettes. Two contours 7 and 

6 arc incompatible if they share some (d — 2)-dimensional face. In this case wc denote 79^^. 

Denote by the set of all contours in the volume A. The subset C N^'^ of compatible 
configurations is defined as 

A"^ = {7? G {0, l}^'' ; 77(7) = if 7 7^ n (2-1) 

that is, a configuration of contours is compatible if it does not contain two incompatible contours. 
We denote X = X^'^ and G = G^'^ . As usual, we endowed X with the product topology. 

For each fixed P G M+, a parameter usually called the inverse temperature, and for each 
finite A define the measure /i^ on by 

exp(-/3X;w-y)=i ItI) 
M^W = ^ ^ ^ (2.2) 

where I7I is the number of plaquettes in 7 and is a renormalization constant making /x^ a 
probability. 

2.2 The loss network 

We introduce a Markov process called (fixed routing) loss network in the set of compatible 
contours. This process was introduced by Erlang in 1917 (see Brockmeyer, Halstrom and Jensen 
(1948), p. 139). An account of its properties can be found in Kelly (1991). In the traditional 
interpretation a contour 7 represents the route taken up by a call. The plaquettes encompasing 

7 are the circuits held by the call. For a finite or infinite set A C Z*^ and / real continuous 
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function on X^, the generator of the process is defined by: 

A'^m = E ^'"^^^ + G -^^1 [fiv + - fiv)] + E ^(^) [/(^ - -^7) - fiv)] (2.3) 

where S.y{6) = 1{9 = 7} and the sum of configurations is defined pointwisely in N*^"^: (f?+'^)(7) = 
77(7) + ^(7)- 111 words, each contour 7 attempts to appear at rate e~'^l7l i^^t it does so only 
if it is compatible with all present contours. Present contours disappear at rate 1. Our first 
result is the following sufficient condition for the existence of a process with generator for 
any (infinite) A. 

Theorem 2.4 // 

a(/3)=sup^ V 10|e-^l^l<oo (2.5) 

then for any ( infinite ) A the Markov process with generator exists and admits at least one 
invariant measure. 

This theorem is proven in Section ^. We denote r;^'^ the corresponding process in A with initial 
configuration C,. We omit the volume superindex when A = 7/: 

4 = ^f^' A = A^' (2.6) 

2.3 Results on the invariant measure 

We say that / has support in T C Z'^ if / depends only on contours intersecting T (not 
necessarily contained in T). Let |Supp(/)| = min{|T| : / has support in T}. When we write 
Supp(/) we mean some T such that |T| = |Supp(/)| and / has support in T. For instance, if 
/(?]) = 77(7), Supp (/) may be set as {x} for any € 7. The results below work for any such 
choice, so one can take the most favorable one in each case. Let 

13* solution of a(/3) = 1, (2.7) 
«o(/3) = El7|e-^l^l- (2.8) 

730 
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Let be some fixed norm, for instance the one given by the Manhattan distance (|x| = 
Yli=i sum of the coordinate lengths). Let the corresponding distance between two 

subsets of Z'^ be 

d(A,T) = min{|x -y| : X G A,y G T}. (2.9) 

Theorem 2.10 If (3 > (3* (that is, a{(3) <1) then the following statements hold: 

1. Uniqueness. For any A C Z'^ there is a unique process rj^^ with generator . The process 
has a unique invariant measure denoted fi^ . For A finite this measure is precisely {2.i). 
For A = Z"^ we denote rjt = rff^, /U = ^u^''. 

2. Exponential time convergence. For any A C and for measurable f on X^, 



sup 



//-E/(r?f'«) < 2||/|U|Supp(/)|^e-^* (2.11) 



where p = i— ^ . 

' 2— a 



3. Exponential space convergence. Let A be a (finite or infinite) subset of 7/ and f a mea- 
surable function depending on contours contained in A. Then 

< 2||/|UaoM2 e-*'^3'^(W,A^) (2.12) 

xeSupp (/) 

where Ms = (1 - a(/3))-i and M3 = (/3 - ]5), for any /3 G (/?*, /3). 

4- Exponential mixing. For measurable functions f and g depending on contours contained 
in an arbitrary set A C Z'^.- 

l/x^(/5)-//Al < 2||/|U||ff||oo(M2)2 Yl |x-y|e-*^3|x-.l (2.13) 

xeSupp {/), 
j/eSupp (3) 



where M2 ad M3 are the same of l \2.1^ ) 



5. Central limit theorem. Let f be a measurable function on X with finite support such that 
pLf = and /^d/P^'^) < 00 for some (5 > 0. Let Tx be the translation by x and assume 
D = Ylx f^if'^xf) > 0. Then D < 00 and 

1 



(0,") p-'^) 



where the double arrow means convergence in distribution. 



This theorem strengthens the results R1-R5 of Fernandez, Ferrari and Garcia (1998). In 
that paper /?* was replaced by a value Pm defined as the solution of qo(/3a/) = ^/{d — 1)- This 
value is strictly bigger than Item 5 generalizes (the central limit) Theorem 7.4 of Dobrushin 
(1996a), where only functions depending on a finite number of contours are considered. The 
above theorem will be proven in Section ^. 

Finally, we prove a Poisson approximation. Consider the equivalence relation induced by 
the translation of contours. Let Gj be a the set formed by one representative containing the 
origin from each class of equivalence of contours with length j. For each Borel set V C M"^ and 
a G M let 

-x-l/2 x + 1/2- 



V -0 = Ixel'^ 



CV). (2.15) 
a a i J 

Fix a contour length j > 0. For each 7 E Gj, let M^^p be defined by 

M,Ay)= E ^^(^-^) (2.16) 

where r/^ is distributed according to the invariant measure for the given (5. Let (My^oo) 7 £ Gj) 
be a family of independent unit Poisson processes in W^. 

Theorem 2.17 For each contour length j, it is possible to jointly construct the d- dimensional 
processes {M^^^ '■ 13* < (3 < 00, 7 E Gj} in such a way that for all regions V equal to a product 
of intervals, 

¥{m^AV)^M^^UV)) < c(|7|,F)(Q(/?) + e-^H/'^) ~ exp(-/3min{2d,|7|M). (2.18) 

where c{\j\,V) is a computable constant. As a consequence {M-y^^ : 7 E Gj} converges in 
distribution to a family consisting of \ Gj \ independent Poisson processes with mean e~^^ . 



Theorem 2.17 is proven in Section |8|. 
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3 Graphical representation of loss networks 

We construct the loss network as a function of stationary marked Poisson processes (a la Harris) 
each of which indicates the attempted birth-times of a contour. A lifetime is associated to each 
attempted birth. The triple (contour , attempted birth , lifetime) is called a cylinder. The loss 
network is constructed by erasing cylinders which at birth violate the exclusion condition. The 
crucial point in this construction is the association of the lifetime to the birth-time. This allows 
us to study the process backwards in time by studying a Markovian oriented percolation process 
of cylinders. In contrast, the standard construction uses independent Poisson processes for the 
birth-times and death-times respectively. In this case the backwards construction looks hard. 

3.1 Marked Poisson processes 

To each contour 7 G G we associate an independent (of everything) marked Poisson process 
with rate e~^\'^\. We call Tkij), 7 G G, the ordered time-events of Nj with the convention that 
To (7) < < ri(7). For each occurrence time 1^(7) of the process we choose an independent 
mark Si{^) exponentially distributed with mean 1. At the Poisson time-event Ti^-y) a contour 
7 appears and it lasts Si{j) time units. 

The random family C = {{(7, 1^(7), /S^ (7)) : i G Z} : 7 G G} consists of independent marked 
Poisson processes. A marked point {'y,Tk{'j), Sk{'j)) G C is identified with 7 x [Tfe (7) , (7) -|- 
<S'fc(7)], the cylinder with basis 7, birth-time Tk{'y) and lifetime Sk{'y)- The life of the cylinder is 
the time interval [Tj^ (7) , Tj^ (7) -|- Sk{'y)]. For a generic cylinder C = (7, t, s), we use the notation 

Basis(C) = 7, Birth (C)=t, Death (C) = t + s , Life (C) = [t,t + s]. (3.1) 

We define incompatibility between cylinders C and C by 

C T^C if and only if Basis (C) / Basis (C) and Life (C) n Life (C) / 0, (3.2) 

otherwise C ~ (7 (compatible). We say that two sets of cylinders A and A' are incompatible if 
there is a cylinder in A incompatible with a cylinder in A': 

A 7^ A' if and only if C / C" for some C G A and C G A'. (3.3) 
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Let S := {Sf{9) : 6 £ G,i > 1) he a (countable) family of iid exponential times of mean 
1 independent of C. These are the lifetimes that, when necessary, will be associated to the 
contours of the initial configuration. Indeed, we identify S with the set of cylinders {{9, 0, S^{9)) : 
9eG,i>l}. For ^ € N<^, let 

s(e)= U [JiiOAs^m} (3.4) 

egG i=l 

the family of cylinders associated to the initial configuration ^, all with birth-time 0. Notice 
that ^ may have more than one cylinder with the same basis. For s < t define 

C[s, t]:={C eC: Birth (C) G [s, t]}. (3.5) 

the set of cylinders born in the interval [s,t]. 

Remark. In this paper we will work with the probability space given by the product of the 
spaces generated by C and S. We call it {Q,J^,F). We write E for the respective expectation. 
In section 4.5.2| we use the direct product of this space by itself, while in Section |5| we need 



to consider countable products of this space. We use the same notation P and E for the 
corresponding probability and expectation in these enlarged spaces. 

3.2 The free network 

For ^ G N*^ define 

Cf(7)= Yl l{Basis(C)=7,Life(C)9i}. (3.6) 

CGC[0,i]US(C) 

The above process, called the free network, is a product of independent birth-and-death processes 
on with initial configuration ^ whose generator is given by 

A'fiO = E ^"'""'[/(^ + '^7) - fiO] + E ^(^)[/(^ - '^7) - fiO]- (3.7) 

The invariant (and reversible) measure for this process is the product measure fi^ on N*^ with 
Poisson marginals 

/{e(7) = k} = ^^^^ exp (e-^l^l) . (3.8) 
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In terms of loss networks, is the process for which all the calls are accepted, that is (7) 
is the number of calls on route 7 at time t when there is no restriction on the number of calls 
a circuit can accept and the initial configuration of calls is ^. 

3.3 Finite- volume construction of a loss network 

In the construction of a loss network in a finite volume A with an initial condition ^ G we 
use only the finite set of Poisson processes (A'^^ : 7 C A) and the finite family of initial lifetimes 
{S^{9) : e C A). Let = {C G C : Basis (C) C A} and S(C), defined as in Q, be such 
that all its cylinders are mutually compatible. We realize the dynamics r/^'^ as a (deterministic) 
function of and S((^). 

We construct inductively K^[0,t], the set of kept cylinders at time t. The complementary 
set corresponds to erased cylinders. At time zero we include all cylinders of S(C) in K^[0,t]. 
Then, we move forward in time and consider the first Poisson mark: The corresponding cylinder 
is erased if it is incompatible with any of the cylinders already in K^[0, t], otherwise it is kept. 
This procedure is successively performed mark by mark until all cylinders born before t are 
considered. Define r/^''' € as 

7]^''^{j) = l{Basis (C7) = 7, Life (C) 9 t} (3.9) 

CGK^[0,t] 

= l|7 G {Basis (C) : C € Kj[0,t], Life (C) 9 t}}, (3.10) 
that is, r/^'^ signals all contours which are basis of a kept cylinder that is alive at time t. We 



show in Section 3A that r/^'^ has generator defined as in (|2.3| ) restricting the sums to the 
set of contours contained in A. It is immediate that fx^ defined in (|2.2| ) is reversible for this 



process. Since we are dealing with an irreducible Markov process in a finite state space, r/^''' 
converges in distribution to for any initial configuration This in particular implies that 
fi^ is the unique invariant measure for this process. Later in the paper we determine the speed 
of convergence. 

Using the same C and S in the construction of r/^'^ and , we have that if ^(7) < ^(7) for 
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all 7 C A, then 

vt'Hl) < Ct^(7), for all 7 C A, (3.11) 

because in the free network all cylinders are kept. 

Since A is finite, there exists a sequence of random times ti = tj(C^) with ti ±oo as 
i zboo such that ^tiil) = for all 7 € A. We can, in particular, consider ti as the entrance 
times of in the set : .^(7) = 0, for all 7 G A}. Since this process has a unique invariant 
measure which gives positive probability to this set, (ti) is a stationary renewal process with 
inter-renewal time with finite mean. We extend the construction of a set of kept cylinders to 
t E M forgetting the set S(C), by doing the above procedure in each time interval with 
the cylinders of C[ti,tj+i]. This can be done because no cylinder intersects {tj : i G Z}. Let 
us denote the resulting set of kept cylinders and r/f^ its projection in the sense of (p.S|). By 
construction has a time-translation invariant distribution. The process rj^ has generator 
and distribution independent oft, hence given by /i^. This implies that for any / : {0, 1}*^"^ — > M 
and any t G M, 

// = E/(r?,^). (3.12) 

Since ??^(7) < ^til) for all 7 G A and ^4(7) has Poisson distribution with mean e~^^"'^, we have, 
taking /(r?) = 1/(7) in ( p. 12]) , that 



/{t? : 77(7) = 1} < e-'^H. (3.13) 

3.4 Infinite-volume construction. Backwards oriented percola- 
tion 

If we try to perform an analogous construction in infinite volume we are confronted with the 
problem that there is no first mark. To overcome this we follow the original approach of Harris 
(1972) [see also Durrett (1995)] and introduce the notion of percolation. The goal is to partition 
the set of cylinders in finite subsets to which the previous mark-by-mark construction can be 
applied. 
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We come back to the infinite- volume construction of Section 3.1. For an arbitrary space-time 



point (x, t) define the set of cyhnders containing the point (x, t) by 

Ai'* = {C € C ; Basis (C) 9 x, Life (C) B t} . (3.14) 

For any cyhnder C define the set of ancestors of C as the set of cyhnders born before C that 
are incompatible with C: 

Af = {C' eC;C' ^ C, Birth (C") < Birth (C)} (3.15) 

xGBasis (C) 

The definition of ancestor of C does not depend on the lifetime of C. Recursively for n > 2, the 
nth generation of ancestors of (x, t) is defined as 

A^'* = {C" : C" G Af' for some C £ A^lj, (3.17) 

and for a given cylinder C, 

A^ = {C" : C" G Af for some C G A^_i}. (3.18) 

We say that there is backward oriented percolation in C if there exists a space-time point 
{x, t) such that A^' * ^ for all n, that is, there exists a point with infinitely many generations 
of ancestors. Let the clan of the space-time point (x, t) be the union of its ancestors: 



n>l 



U A^'- (3-19) 



In the next theorem we give a sufficient condition for the existence of the infinite-volume 
process in any finite time interval in terms of backwards percolation. 

Theorem 3.20 // with probability one A^'* n C[0,t] is finite for any x G Z'^ and t > 0, then 
for any (possibly infinite) A C Z"^, the process with generator is well defined for any initial 
configuration ( G X^and has at least one invariant measure /x^. 

Proof. We construct the process for h. = TL'^\ the construction for other regions A is analogous. 
The initial configuration is denoted C, ^ X and the initial cylinders are given by S(C), defined 
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in ( [j.4D . Note that that all the cyhnders of S(C) are mutually compatible. We then partition 
S(C) U C[0,t] into a set of kept cylinders, denoted by K, and a set of erased cylinders denoted 
by D. 

The construction is as follows. First, all cylinders in S(C) are kept. Second, for each a; E A 
the percolation clan of {x,t) in [0,t], A^'* n C[0,t], is partitioned in kept and deleted cylinders 
as in the finite- volume case. To do so we order the cylinders of A^''* n C[0,t] by birth-time. 
This can be done because by hypothesis A^'* n C[0, t] has a finite number of cylinders. Then we 
successively classify each cylinder as kept if it is compatible with all cylinders already classified 
as kept (including those in S(C)); if not, we classify it as erased. We denote the resulting sets 
K^'*[0, t] and D^'*[0, t], respectively. 

Denoting 

K^O,*] := |jK^'*[0,t], Df[0,t] := [Jnfm, (3.21) 

we have that 

K(_-[0,t]UD^[0,t] = C[0,t] U S(C) . (3.22) 

Indeed, the classification of any given cylinder C G C[0, t] depends only on (a) its ancestors in 
[0, t], A*^ n C[0, t], and (b) on the finite subset of S(C) of cylinders which are incompatible with 
some of the ancestors of C in [0, t] . Therefore there is no inconsistency: K^'* [0, t] n D^'* [0,t] = 
for all X / y and K^'*[0,t] C K^''*'[0,t'] for t < t', if {x,t) G C for some C G K^''*'[0, t']. 
The process is now defined as in (|3.9D by 

7,^(7) = e {Basis (C) : C G KjO,t], Life (C) 3 t}}. (3.23) 

The reader can check that for finite A the above construction is equivalent to that of Section 
|3.3| . Applied to the set of cylinders of C^[0, t] it yields the set K^[0, t] defined in the paragraph 
preceding formula ( |3.9D . 

To show that r]^ has generator A, denote rjt = rfi and K = and write 

[fim+h) - fivt)] 

l{Birth(C) G [t,t + /i]}[/(r/t + 5Basis(C)) -/(^t)] 

C&K[0,t+h] 
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+ Yl liLife {C) 3 t, Life (C) ^ t + h}[f{rjt - ^Basis (C)] 

CGK[0,t] 

+ {other things}, 



fim)] 



(3.24) 



where {other things} refer to events with more than one Poisson mark in the time interval 
[t, t + h] for the contom's in the (finite) support of /. Since the total rate of the Poisson marks 
in this set is finite, the event {other things} has a probability of order /i^. Now, denoting 



N^it,s) = #{k:Tkij)eit,s)} 



(3.25) 



we have 



^ l{Birth(C7) G + h]} 1{C e K[0,t + h]} [/(r?* + 5Basis(C)) " fiVt)] 



cec 



Yl l{^Basis (C) [t, t + h] = l} l{Basis (C) ~ Basis (C), VC" G K[0, t] : Life (C) 3 t} 
Cec 

x[/(^i + <^Basis(C)) - fiVt)] 



Y n^^t, t + h] = l} l{vt +5^eX} [fiijt + 5^) - fiTjt)] 



(3.26) 



To compute the second term of ( 3.24 ), observe that Life (C) is independent of Birth (C) and 
both the event {C £ K[0,t]} and rjt are .Ft-measurable. Here Tt is the cr-algebra generated by 
the births and deaths occurred before t. Hence 



PILife (C) 9 t, Life (C) ^ t + h 



Life {C)^t^h Life {C) 3 t] (3.27) 



and 



E 



1{C G K[0,t]} l{Life (C) 9 t, Life (C) ^ t + /i} - 5Basis(C)) " /(??*)] 

c 



E 



c 



[^P(^Life(C) + Life((:7) 9 1{C G K[0,t], Life (C) 9 t} 

x[/(??t-5Ba.i.(C7))-/(^t)]l- (3.28) 



Since Life (C) is exponentially distributed with mean 1, 



»(^Life(C) ^ 



Life (C7) 9 t 



h + o(/i) 



(3.29) 
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Taking the expectation of (|3.24|) and substituting ( 3.26 )-( 3.29 ) we get 



= Y.^ e-'^l^l ^{l{r,t + 5^(^X} [fir^t + S^) - fivt)]) + o{h) 

7 

+ Y^hEUh)[fivt - S^) - fivt)]) +oih) (3.30) 



7 

which dividing by h and taking Umit gives 

^ = Am4)- (3.31) 

The existence of an invariant measure follows by compactness as our process is defined in 
the compact space X. See Chapter 1 of Liggett (1985). □ 

We show in the next theorem that under stronger hypothesis the process can be constructed 
for times in the whole real line. Since the construction is time-translation invariant, the distri- 
bution of r]t will be invariant. 

Theorem 3.32 // with probability one there is no backwards oriented percolation in C, then 
the process with generator A can be constructed in ( — cxd, oo) in such a way that the marginal 
distribution of rjt is invariant. 

Proof. The lack of percolation allows us to construct a set K C C as K^[0, t] was constructed 
from C[0,t] U S(C) in the proof of the previous theorem. We just proceed clan by clan and 
simply ignore the cylinders of S. Note that K is both space and time translation-invariant by 
construction. Analogously to the previous theorem we define r]t as the section of K at time t: 



Vt 



(7) = ij^ e {Basis {C):C£ K, Life (C) 3 t}}. (3.33) 



By construction, the distribution oir]t does not depend on t, hence its distribution is an invariant 
measure for the process. □ 

Let us denote fi the distribution of r]t, in anticipation of the fact that this is precisely the 



measure of Theorem 2.10 
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As in the finite case, 

Vtil) < 6(7) (3.34) 
for all 7 € G. This implies that the distribution /x inherits property ( p.l3| ): 

KvH7) = '^) = E^t(7) < mh) = e-'^l^l. (3.35) 

Let 

A(T) = U^exA^'O (3.36) 

be the clan of T C Z'^ (at time 0). 

Remarks. 1. It follows from ( |3.33| ) that for any t G M and continuous /, 

fif = Ef{r,t). (3.37) 

2. The presence/absence of a contour 7 at time t depends only on the clan of ancestors of 
(x, t) for any x € 7 through a certain function. More generally, for each / there exists a function 
^> such that 

/(r?t) = cI>(u,.eSupp(/)A^'*). (3.38) 

For instance, r]f(-y) = 1 if and only if A^'* contains a cylinder in K with basis 7 whose life 
contains t. This depends only on the set A^'*. In particular, with the notation ( [3.36 ), 



/(r/o) = $(A(Supp (/))). (3.39) 
Analogous statements are true for the process starting with a fixed configuration at time zero: 

/(r?^) = cI>(A^(Supp(/),[0,t])) (3.40) 

where 

Ac(T,[0,t]) = [(U,6tA"'*) n C[0,t]] U [Ce S(C) : {C} ^ (U,6tA"'*) n C[0,t]} (3.41) 
is the set of cylinders in C[0, t] US(C) which determines the value of /(jyf ) when / has support T. 
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4 Percolation, space-time convergence and mixing 

In this section we exploit the relation between the loss-network process and the absence of per- 



colation to prove a more precise version of Theorem 2.10, In the proof of the mixing properties 



we shall need a continuous-time construction of the backwards percolation clan. 

4.1 The key theorem 

The precise statement of next theorem requires the notion of non-oriented percolation in a time 
interval. For any time interval {s,t) and any space-time point {x,t') define 

G^'*' [s, t] = |C7 G C[s, t] : Basis (C) B x, Life (C) 3 t'} (4.1) 

and 

G^'*'[s,t] = jc G C[s,t] : Basis (C) 9^ Basis (C), for some C G G^'l'js,t]}. (4.2) 

Notice that in the definition of Gn there is no exigency that the birth time of C be previous 
to the birth-time of C or that the lifetimes intersect. Let 

G^'''[s,t] = [jGf[s,t]. (4.3) 

fc>0 

We say that there is no (non-oriented) percolation in [s,t] if for any space-time point {x,t'), 
G^'*'[s,t] contains a finite number of cylinders. We will show later that the condition a < 00 
is sufficient for the existence of an h such that the probability that there is no non-oriented 
percolation in [0, h] is one. 

In addition we need the following definitions. 

• The time-length and the space-width of the family of cylinders A^'* are respectively 

TL (A^'*) = t - sup{s : Life (C) 9 s, for some C G A^'*}, (4.4) 
SW (A^-*) = I U Basis (C) . (4.5) 

In words, the space- witdth is the number of sites occuppied by the projection of the bases 
of the cylinders in the family. The time-lenght is the lenght of the time interval between 
t and the first birth in the family of ancestors of {x,t). 
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• A^(Supp (/)) is the set of ancestors of Supp (/) constructed from as A(Supp (/)) was 
constructed from C. [Notice that this is not the same as A(Supp (/)) fl C^]. 

In item 4 of next theorem we enlarge our probabihty space to the direct product of our 

working space with itself: (fi,^, P) x (Q,J^,¥). As warned before, we continue to use P and E 
for the probability and expectation of this space. 

Theorem 4.6 Assume that there is no backwards oriented percolation with probability one. 

Then, 

1. Uniqueness. The measure ji is the unique invariant measure for the process r)t. 

2. Time convergence. For any function f with finite support, 



lim sup 



lE/(??t^)-M/ =0. (4.7) 



Furthermore, 

/x/-E/(,7^)^ 

< 2||/|UP(u,eSupp(/){A"'VS(C) orTL(A-'*)>t}) (4.8) 



sup 

vex 



< 2||/||oo Yj [p(tL(A^'°) > 6t) +e-M*E(sW(A^'0))] (4.9) 

a;GSupp (/) 

for any b € (0, 1). 

3. Space convergence. As A ^ Z*^, /j,^ converges weakly to ji. More precisely, if f is a 
function depending on contours contained in a finite set A, then 

< 2||/||ooP(A(Supp(/))7^A^(Supp(/))). (4.10) 

4. Mixing. If in addition there exists a value h such that there is no (non- oriented) percola- 
tion in (0, h) with probability one, then for f and g with finite support, 

lim HfT^g) - fif fig\ = (4.11) 

\x\—>-oo 

where Tx is translation by x. More precisely, 

Hf9)-l^fl^9\ < 2||/||oo||5||ooP(A(Supp(/)) / A(Supp(5))) (4.12) 
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where A(Supp(5)) has the same distribution as A(Supp(g)) but is independent of 
A(Supp(/)). 



Existence of /i has been proven in Theorem 3.32. In the rest of the section we prove the 
other properties. 



4.2 Time convergence and uniqueness 

We use the same Poisson marks to construct simultaneously the stationary process rjt and a 
process starting at time zero with an arbitrary initial configuration The second process is 
denoted rj^ (as before). This is what in the literature is known as coupling. By construction 
[cf. (|3:2lD and (pSj)], the process rj^ ignores the cylinders in C with birth-times less than 
but takes into account the set of cylinders with basis given by the contours of the initial 
configuration and birth-time zero, S(C) = {{d, 0, S^{9)) € S : rj{9) = 1}. Recall that the times 
S'^{9) are exponentially distributed with mean 1 and independent of everything. 
By ( lO^ ) and ( p^ ), 



sup 



^f{v^)-^^f 



sup 



(4.13) 



Since we are using C to construct rjt and C[0,t] U S(C) to construct rj^, it follows from ( 3.38 ) 
and ( CT) that 

fiv^)-f{Vt)\ = |$(A^(Supp(/),[0,t])) -<I>(u,eSupp{/)A^'*^ 

< 2||/|Ul{(s(C)9^U,6Supp(/)A^'*) orTL(A-'*)>t} (4.14) 

To see this notice that {A^'* C C[0,t]} = {TL (A^'*) < t}, and that U^gSupp (/) A^'* ~ S(C) and 
U.gsupp{/)A^'* C C[0,t] if andonly if U,esupp{/)A^'* = A^(Supp (/), [0, t]) . Display ^ 



shows (4 



To prove the weak convergence ( ^.Tj ) we fix 6 G [0, 1] and bound the indicator function in 
the right hand side of (4.14) by 



l{tL(A^'*) > btj + i|tL(A^'*) < bt,A^'^ r/. S(C)}. 



(4.15) 
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The expected value of the first term in (|4.15| ) goes to zero because A^'* has a finite number of 
cylinders with probability one. The second term in ( [4 .151 ) is bounded above by 

l|max{5?(7) : ((7) = 1 and 7 Basis (C) for some C £ A^-*} > (1 - (4.16) 

Since S and A^'* are independent and are iid exponentially distributed random variables of 
mean 1, 



x,t 



E(^max{5°(7) : C(7) = 1 and 7 / Basis {C) for some C G A^'*} > (1 - b)t 

= 1 _ (1 _ g-(i-f')*)l{7:C(7)=i and 77^Basis(C) for some CgA^'*}| 

Since ^ is a configuration of compatible contours, it contains at most one contour per site, i.e. 
\{jB X : C(7) = 1}| < 1 for all x G Z'^. This implies that at most SW(A^'*) cylinders of C{C) 
can be incompatible with cylinders in A^'*. Hence, ( 4.17| ) is bounded by 



l_(l_e-(i-'')*)SW(A-0. (4.18) 
The expectation of ( |4.1^ ) is given by 

[1 - (1 - e-(i-^)*)"] P(SW (A^'°) = nj (4.19) 

n>l 

because the distribution of A^'* does not depend on t. Our hypothesis of no backwards 
oriented percolation implies that A^'^ contains a finite number of (finite) contours. Hence 
'}2n>i ^"(SW (A^''') = n) = 1 and by dominated convergence ( 4.19| ) goes to zero as t ^ 00. This 
proves (|4.7|). 



To prove ( [4. 91 ) we start from the expectation of ( [4.15| ) and use (4.19) to bound the expected 
value of the second term by 



k 



e-a-b)t ^ p(^sW (A-'O) =n)Y.{l- e-^'-'^\ 

n>l k=0 

< e-(i-'')*^nP(^SW(A^'°) =n) 

n>l 

< e-(^-'')* EfsW(A^'°)). (4.20) 



The above arguments prove that the process converges, uniformly in the initial configuration, 
to the invariant measure /x. An immediate consequence is that /x is the unique invariant measure. 



This concludes the proof of (1) and (2) of Theorem 4^. □ 
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4.3 Finite- volume effects 



To prove inequality ( 4.101) we use (l3l3^ ), and (|3^ to get 



E 



$(A(Supp(/))j -$(A^(Supp(/))jJ, (4.21) 

where ^ is the function referred to in ( p. 39 ). By definition, 

<I>(a(Supp(/))) < ll/lloo. (4.22) 

Hence, inequahty (4.10) follows from ( [4.21| ). 

Since the spatial projections of the set of ancestors of Supp (/) are finite, the right hand side 
of ( 4.10| ) goes to zero, proving, in particular, the weak convergence of /i^ to /i. □ 



4.4 Mixing. Its relation with a coupling construction 



The proof of (4) of Theorem 4.6 is very similar in spirit to the above proof but it requires a 
somewhat more delicate argument based on the coupling of two continuous-time versions of the 



backwards percolation process. We first notice that (4.11) is a straightforward consequence of 
( [4.12 ), because in the absence of backwards percolation, the spatial projections of the set of 
ancestors of Supp(/) and Supp ((7) are finite. This implies that the right hand side of ( 4.12 ) 
goes to zero. 

To prove ( |4T2D we use ( ^ ) and (|]38|) to get 

= E{f{r,o)g{rio))-Ef{r,o)Eg{r,o) 



E 



$ A(Supp(/))U A(Supp(<7)) 



E 



«l> A(Supp (/)) «!> A(Supp(ff)) 



(4.23) 



where $ is the function referred to in (|3.38D and (A(Supp (/)), A(Supp ((/))) has the same 
marginal distributions as (A(Supp (/)), A(Supp (5))) but its marginals are independent. 

Identity ( [4.23| ) shows that to obtain ( [4.11| ) it is enought to construct a coupling (joint 
construction) of the four processes 

A(Supp (/)) , A(Supp (g)) , A(Supp (/)) , A(Supp (g))) , 
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such that 



A(Supp(/)) = A(Supp(/)) 



(4.24) 



and 

A(Supp(5)) ~ A(Supp(/)) impHes A(Supp (5)) = A(Supp (5)) • (4.25) 

Indeed, from ( 4.22| ) and ( 4.2^ )-( |4.25| ) we obtain that the last hne of (4.23) is bounded above by 
the right-hand side of ( 4.121 ). 

In the remaining of the section we discuss the construction of the couphng with properties 
( [4.24D -( [l.25| ). The construction is natural and straightforward, but unavoidably technical. As 
an alternative we mention the approach based on "duplicated variables" (von Dreifus, Klein 
and Perez, 1995; Bricmont and Kupiainen, 1996), which is probabilistically simpler but requires 
some combinatorial input. 



4.5 Construction of a four-clan coupling 

We need to couple two clans in the same random set of cylinders with two independent copies 
with the same marginal distribution. Moreover, to strengthen our results we need to ensure that 
the marginal realizations remain the same as much as possible. The coupling (Section |4.5.2 ), 
is based on a construction of backwards percolation clans as non-homogeneous continuous time 
Markov processes (Section fl.5.1| ). The hypothesis on the absence of non-oriented percolation 
for some time interval (0, h) is needed for the infinite-volume construction of the coupling. 

4.5.1 A continuous-time construction of the backwards percolation clan 

For T CZ'^ define 

At(T) = |C7' G A(T) : > Birth (C) > = A(T) n C[-t, 0], (4.26) 

that is, the set of cylinders in ^(T) with birth-time posterior to —t. The inclusion of a new 
cylinder in the time interval [t, t + h] depends on the existence of a birth Poisson mark in 
[—t — h, —t] whose corresponding cylinder is incompatible with some C € Aj(T). That is, if C 
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is a cylinder with Basis (C) ^ Basis (C) for some C" G A^, 



Kt^h = A U C 



At = A, At/ = At/, t' G [0,t) 



= P{C G C : Birth (C7) G [-t - /i, -t\ , Death (C) > t - TI (A, Basis (C))} + o(/i) . 

(4.27) 

We have denoted 

TI (At, 7) = min{Birth {C') : C G At, Basis {C) / 7} (4.28) 

and abbreviated At(T) = At. The remainder o{h) is the correction related to the probability 
that C is not the only cylinder born in [—t — h, —t]. Since the birth-time is independent of the 
lifetime which is exponentially distributed with rate one, 

At+h = A U C I At = A, At/ = At/, t' G [0,t) 
= P|C7 G C : Birth (C) G [-t - h, -t]} P(^Life (C) > t - TI (A, Basis (C))) + o{h) 

^ ;^g-/3|Ba.i.(C)|g-t+TI(A,Ba.i.(C)) ^ ^ (4 29) 

This implies that when the configuration at time t~ is A, a new cylinder with basis 7 is included 
in At(T) at rate 

From (|4.2S| ), as in the computation of the forward Kolmogorov equations, we get 
'4f(At) 

It-Tl (At ,7) 

where the sum is over the set {7 G G : 7 9<^ Basis (C) for some C G At}. This equation 
characterizes the law of the process At(T) as a non-homogeneous Markov process. 

We now construct At(T) by combing the Poisson marks backwards in time in a continuous 
manner. 

Finite- volume case. If we only consider contours contained in a finite set A C Z*^, there is 
only a finite set of possible bases for the cylinders and the Poisson marks are well ordered with 
probability one. The construction proceeds mark by mark backwards in time. Set Ao(T) = T. 
If there is a Poisson birth mark at time —t whose corresponding cylinder is called C", then 
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E 



A„ < s < t ) = V / ds e-' e-^l^l /(At U (7, t, s)) - /(At) (4.31) 



• if C" / C for some C G Ai_(T), set Ai(T) = At_(T) U {C"}. 

• if C" ~ C for all C € Ai_(T), set At(T) = Aj_(T). 
where the incompatibility between cylinders was defined in (|3.2D. 



Infinite volume case. In infinite volume the construction can be performed using a perco- 



lation argument as in Section 3.4. By hypothesis, there exists an h such that each cylinder born 
in the interval [—h, 0] belongs to a finite non-oriented clan. Hence the set of cylinders born in 
the interval [— /i, 0] can be partitioned in connected families: 

C[-h,0]=]JUk[-h,0] (4.32) 

fc>0 

where the sets Hfc[s,t] are the maximal sets of cylinders with the property that cylinders in 
different H^'s are compatible. We can then well-order the birth-time of the cylinders inside each 
Hfc and proceed as for the finite- volume case. This yields the process At(T) for t € [0, h]. To 
extend the construction for arbitrary t > 0, we simply repeat the previous procedure in [h,2h], 
[2h,3h], etc. 

4.5.2 A coupling between two interacting and two independent clans 

We take two independent marked Poisson process whose marks and cylinders we respectively 
call blue and red. We enlarge our probability space and continue using P and E for the 
probability and expectation with respect to the space generated by the product of the blue 
and red Poisson processes. Using these marks we construct simultaneously the processes 
(At(Ai), Aj(A2), At(Ai), At(A2)), for Ai,A2 C Z'^, in the following way. 

1. The processes Aj(Ai) and At(A2) are constructed using only the blue marks, as described 
in Subsection 4.5.1, and ignoring the red marks. Hence, they are the clans of Ai and A2 



respectively. 

2. The process Aj(Ai) is also constructed only with the blue marks, hence it coincides with 
At(Ai). 
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3. The process At(A2) is constructed with a precise combination of blue and red marks in 
such a way that (a) it coincides with At(A2) for a time interval that is as long as possible; 
(b) it is independent of Af(Ai), and (c) it has the same marginal distribution as At(A2). 

Property (^) is achieved in the following way. 

Finite- volume case. If both Ax and A2 are finite sets, we order the marks by appearance 
and introduce a flag variable, Flag(t) G {0,1}, which indicates if some cylinder of At(Ai) is 
incompatible with some cylinder of At(A2): 

Flag(i) = C" for some C G Ai(Ai), C" G At(A2)}. (4.33) 

We now proceed as follows, mark by mark backwards in time. First, we set Flag(O) = 0. The 
construction guarantees that Flag(t) = implies At(Aj) = At(Aj) for i = 1,2. 

• If at time —t a (blue or red) mark is present and Flag(t— ) = 0, then 

— If the mark is blue and the corresponding cylinder can be included in A(_(Ai) but 
not in At(A2), then include it in At(Ai) and A((Ai). Analogously, if it can be 
included in At_(A2) but not in A((Ai), then include it in A((A2) and A((A2). Keep 
the Flag = 0. 

— If the mark is blue and the corresponding cylinder can be included in both At„(Ai) 
and A(_(A2), then include it in both At(Ai) and A((A2) but include it only in 
At(Ai). Set the Flag = 1. 

— If the mark is red and the corresponding cylinder can be included in both At_(Ai) 
and Kt~{^2), then include it only in Ai(A2). Set the Flag = 1. 

— If the mark is red and the corresponding cylinder can be included in either At_(Ai) 
or At_(A2) but not in both of them, then ignore the mark. Keep the Flag = 0. 

— If the mark is red or blue but the corresponding cylinder can be included in neither 
A(_(Ai) nor At„(A2) then ignore the mark. Keep the Flag = 0. 
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• If at time —t a (blue or red) mark appears and Flag(t— ) = 1, then use blue marks for 
At(Ai), At(A2), and At(Ai) and red marks for At(A2). Keep the Flag = 1. 

To verify that the above coupling has the right marginals, it suffices to notice that the rate 



of inclusion of cylinders in each one of the marginals is precisely given by (4.30). 



Infinite- volume case. In the case in which at least one of Ai and A2 is infinite, we consider 
families analogous to the given in ( B) but defined using the time interval [—h/2, 0] and 
both red and blue marks. Therefore in the combined set of cylinders there is no non-oriented 
percolation and we can construct the coupling working in a finite set of H/j's at a time. We 
then continue working in time intervals of lenght h/2 to reach arbitrary times t. 

By construction, the coupling satisfies ( 4.24] ). It also satisfies ( [4.25|) because the flag changes 



from to 1 (and remains 1 forever) the first time a cylinder of At(Supp(5)) is incompatible 
with a cylinder of Aj(Supp(/)). This implies 

Flag(oo) = i{a(Supp(/)) ^ A(Supp(5))}. (4.34) 

5 Branching processes. Time length and space width 

In this section we estimate the time-length and space-width of the families of ancestors A^'*. We 
follow the well known approach of introducing a branching process that dominates the backward 
percolation process (see eg. Hall, 1985), though we must consider multitype branching. The main 
result of this section is Theorem ^.27| which shows that the hypotheses of Theorem 2.10 lead to 
exponential upper bounds of both TL (A^'*) and SW (A^'*). 

5.1 Multitype branching processes 

We introduce a multitype branching process B„, in the set of cylinders, which dominates A„. 
To do this we look "backwards in time" and let "ancestors" play the role of "branches" . In par- 
ticular, births in the original marked Poisson process correspond to dissapearance of branches. 
We reserve the words "birth" and "death" for the original forward-time Poisson process. 
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We start by enlarging our probability space and defining, for any given set {Ci, . . . ,Cfc}, 
independent random sets ' with the same marginal distribution as A^^ ' . The important point 
here is that 

k k 

|JAf»c|JBf^ (5.1) 

i=l i=l 

The proof of this fact relies in fixing a way to distribute common ancestors. For example, 
consider the total order ^ in the set of cylinders induced by the birth-times. That is C ^ C" 
if and only if Birth (C) < Birth (C). For any finite set of cylinders {Ci, . . . , (7^} such that 
Ci -< Cj_|_i, i = 1, . . . , fc — 1, define 

A?=A?\(^Ua?)- (5-2) 

^ c 

This ensures that are independent sets and 

k k 

UAf" = UAf'. (5.3) 

i=l i=l 

On the other hand, for any C, Aj^ is stochastically dominated by Af [that is, there exists a 
joint realization (Af , Aj^) such that P(Af C Af ) = 1]. From this observation and ( ^.3| ) we get 



The procedure defined by Bi naturally induces a multitype branching process in the space 
of cylinders. We define the n-th generation of the branching process by 

B^ = {Bf : C G B^_ J (5.4) 

where for all C , Bf" has the same distribution as Af" and are independent random sets 
depending only on C . Inductively, 

A^ C B^. (5.5) 

Indeed, 

U <= U (5.6) 

where in the definition of Aj^' we use {Ci, . . . , C^} = U"jQ^Ap. Hence, the inductive hypothesis 
Af C Bp , for i = 1, . . . , n - 1, yields (^ . 
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In consistence with our previous notation, we denote 

= IJ ; B-'* = IJ B^'* ; B^ = U B^'", (5.7) 

n>0 n>0 xeX 

the branching clans of C, {x,t) and T (at time 0) respectively. By (|5.5|), 



C B^; A^'* C B^'*; A^ C B^. (5.8) 

Defining the time length and space width of this clan as in ( [4.4| ) and ( |4.5| ), we get 

TL (A^) < TL (B^) ; TL (A^-*) < TL (B^-*) ; TL (A^) < TL (B^) , (5.9) 

and similarly for the respective space widths. 

The (multitype) branching process B„ induces naturally a multitype branching process in 

the set of contours. For a cylinder C with basis 7 and birth-time 0, define G N*^ as the 
number of cylinders in the nth generation of ancestors of C with basis 9: 



hl{9) = jC e B^ : Basis (C) = 9^ . (5.10) 



This process will be useful in estimating the space properties of the clans of ancestors. We have 
the following relationship: 

Y,K{9) = \B^\. (5.11) 

e 

The process b„ is a multitype branching process whose offspring distributions are Poisson 
with means 



POO 

m{-f,9) = l{j^9}e-^\^\ e'^dt 

Jo 

= l{7 9<i0}e-^l^l. (5.12) 



To see this, notice that the cylinders C with basis 9 that are potential ancestors of C (with 
basis 7) form a Poisson process of rate 1(7 9*^ 9} e^'^'^L Each of those cylinders is an ancestor of 
C if its lifetime is bigger than the difference between the birth-time of C and C . The lifetimes 
of different cylinders are independent exponentially distributed random variables of rate 1. The 
probability that the lifetime of any given cylinder is bigger than t is given by e~*. Hence, the 
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birth-times of the ancestors of C with basis 6 form a (non homogeneous) Poisson process of rate 
depending on t given by 1{7 / 9}e-^\^\e~K The mean number of births is therefore given by 



Lemma 5.13 The means ( 5.11 ) satisfy 



^m''{j,e) < ^|0|m"(7,0) < |7|a", (5.14) 



where a is defined in (2.5) 
Proof. 



m 

71:71/7 72:72/71 6:6'/7„_i 



I7I I71I ^-^ I7, 

71:71/7 72:72/71 0:6l/7„_i 
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n-1 



< |7|('sup ^ l^e-^l^l") .□ (5.15) 

^ ^ 61:61/7 ^ 

This lemma shows, in particular, that the branching process b„ is subcritical if a < 1. 



5.2 Continuous- time branching process 

Let C be a cylinder with basis 7 and birth-time 0. Combing backwards continuously in time the 
branching clan B*-" we define a continuous-time multitype branching process ipt{0) = number of 
contours of type 9 present at time t (of this process) whose initial configuration is 6^. Each C £ 
B*^ is a branch, that is, belongs to the first generation of ancestors of a unique cylinder U{C') 
in B*-" . In the branching process "0* ah the branches (ancestors) of U{C') appear simultaneously 
at the birth of U{C'), that is when U{C') dissapears if we look backwards in time. Therefore 
the part of C in the interval [Birth (C/(C")), Death (C")] is ignored. Formally, 

^7(0) = G B^ : Basis (C) = 9, Birth (C) < -t < Birth ([/(C)), Life (C) 3 t} . (5.16) 
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In the process V't; each contour 7 hves a mean-one exponential time after which it dies and 
gives birth to kg contours 9, 9 £ G, with probability 



il k;i ('-'^^ 

for kg > 0. These are independent Poisson distributions of mean m{'j,9). The infinitesimal 
generator of the process is given by 

^ _ ^ _ „ pm.(-r,e) ( o\ri{0) 

i^/W = E^W E n lfi^ + ri-S,)-fm (5.18) 

where ip,T] e yo = {ip e N^; J2e ^(^) < ^'^'^ ^0(7) = {V' G 3^0; ipid) > 1 implies 6* / 7} and 
/ : 3^0 ^ N. 

The branching process ip'l allows us to estimate the time-length of a clan, due to the obvious 
fact: 

^jj2{9) = implies TL (B^) < t . (5.19) 

e 

Let Mt{'~f, 9) be the mean number of contours of type 9 in ipt and Rtij) its sum over 9: 

Mt{j, 9) = ¥.^2 {9) ■ Rt{j) = Mt{7, 9) . (5.20) 

e 

The bound we need is given in the next lemma. 

Lemma 5.21 The mean number of branches, Rtij) satisfies 

P(EV'7(^) > 0) < Rtil) < |7|e("-i)*. (5.22) 
e 

Proof. The first inequality is immediate because V'7(^) assumes non- negative integer values 
and Rtil) is its mean value. 

To show the second inequality we first use the generator given by (|1|) to get the Kol- 
mogorov backwards equations for Rt{^)- 

d 



dt 

Since i?o(7') = 1; the solution is 



Rt{l) = Ym{-i,-f')Rt{i) - Rt{l). (5.23) 



exp[t(m — /)] 
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(7,7') (5.24) 



where m is the matrix with entries 771(7,7') ^ is the identity matrix. This can be rewritten 
as 

fn 

n>0 ■ 7' 

n>0 



where the last bound is just the leftmost inequality in (5.14). □ 



5.3 Time length and space width 

We are now ready to provide bounds for the time length and space width of the percolation 
clan. 

Theorem 5.27 If f3 > p* (i.e. a{p) <1), then 

(i) The probability of backward oriented percolation is zero, 
(a) For any positive b, 

P(^TL(A^'*) > < aoe"(^-")''* (5.28) 

(Hi) 

ao(/3) 



E 



(SW,A.^.)) < ^ (5.2.) 

(iv) 



(v) 



E(exp[aSW(A-)]) < ^^^^ (5.30) 



SW(A=''*)>^) < ^^e-(/3-/3)« (5.31) 
1 — a{(3) 



for any (3e (/3*,/3). 
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Proof (i) We follow an idea of Hall (1985). For each C € C we use the domination (5.5) 
and the identity (|5.11| ). Therefore, to prove that there is no backward oriented percolation it is 
enough to prove that, for fixed 7 

P( ^ b;^(0) / for infinitely many n j = 0. (5.32) 

Since b^(0) assumes non negative integer values, by Borel-Cantelli lemma, a sufficient condition 
for ([5^ is 

^^m"(7,0)<cx). (5.33) 

n e 

But this follows from Lemma |5.13| . 

(ii) By ( 5.1S| ) and (|5.9D for each x G 7 and s <t 

^i)J{e) = implies TL(A^'°) < t. (5.34) 
e 

Hence, 

P(TL(A-'0)>t) < Y.nm{l) = ^)Rt{l)< Y.e-^\^\ Rt{i) < aoe-^^"")* (5.35) 



by the rightmost inequality in ( 5.22 ). 

(iii) We find upperbounds for the space diameter of the backwards percolation clan through 
upperbounds for the total number of occupied points by the multitype branching process 



defined by (|5.10D . In fact 



sw(A-'0)< j^7?o(7)EEi^|b;^(^)- 

79a; n 



(5.36) 



By (|]3|), Er/o(7) < e-^l^l, hence by (|1|) 

^(E^oWEEl^l^-W) ^ 5]e-/5H^j;i0|m"(7,e) 



79a; 



n e 



79a; 



(5.37) 



(iv) Write 
E(gasw^ = ^e'^^F(SW =1) 

< E^"'E E l{ l7iU...U7fc|= 41^(71 3 0, b7H72) > 1,..., b^'^-^Tfc) > l) • 

(5.38) 



k 7iv,7* 
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By the Markovian property of b„ we get 

P(7i9 0,b7^(72)>l,...,bf-H7;t)>i; 

= P(7i B 0) nhf (72) >!)••• P(bf - (7fc) > 1) 

k 

< l{ 71 3 0,71 / 72,..., 7fc-i 9^74 n^~^'^'' (^•39) 
Substituting this in ( 5.3^ ) and using that 



1=1 



e'^^l{|7iU...U7fc| < l{|7iU...U7fc| =4 exp(a^|7i|) (5.40) 

i=l 

we get 

k 

E^^aSW) < ^ ^ l{7l^0,7l9^72,...,7fc-l9^7A:}exp(-(/3-a)J^|7^ 

k 7i,.--,7fc «=1 
k 7i3 727^71 7fe7^7fc-l 

< ao{p -a)Y^a{P -a)''. (5.41) 

fc>0 

(v) It suffices to use (iv) and the exponential Chevichev inequahty and to notice that a must 



g-(/3-a)|7fcl 



be less than /3 — /?* to avoid a zero in the denominator of (|5.30| ). □ 

Remarks. 1. Part (ii) can in fact be proven by a more elementary argument not requiring 
the continuous-time construction of Section p^. The argument gives the same rate of decay as 



in ( 5.28 ) but a worse leading constant. Let us sketch it. 



TL (A^'O) >t) <Y, E 9 0, 71 / 72, ... , Ik-i ^ Ik} + --- + Sk>t) (5.42) 

k 7iv.7fe 

where Si are independent mean one exponentially distributed random variables and independent 
of 7j. The time Si represents the period between the birth of 7^ and 7^+1. As the sum of 
independent exponentials is a gamma distribution with parameters k and 1, 

k 

+ • ■ ■ + 5fc > t) = e-* V ^ . (5.43) 

Therefore ( |5.42| ) is bounded by 

<=-' E I E E i^'i ^"^'-'o E 17*^*°' ^ ^ ^ E 

1=0 k>i -yiBO 72/71 7fe/7fc-l 
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2. In Fernandez, Ferrari and Garcia (1998) we offered an alternative proof of part (iii), 
based on a the computation of the exponential moment of the total population of a subcritical 
single-type branching process, which dominates the space width. However this proof works in 
a smaller range of /3. 



6 Proof of Theorem [274 



The following theorem shows that the condition a < oo implies the hypothesis of Theorem 3.20 



This proves Theorem 2.4 



Theorem 6.1 If a < oo, then for all x and positive t the set A^'^[—t, 0] has a finite number of 
cylinders with probability one. 



Proof. Let C be a cylinder with basis 7 and birth-time 0. Recall the definition of U (C) just 
before display (5.16) and define 



V;7(0) = |C7' G : Basis (C) =e,-t< Birth {U{C'))^ . (6.2) 
The process ip] signals all contours born in [0,t] in the process ip^ . Notice that for 2; € 7, 



, x-,Or 



-t,0] 



< 



B^'0[-t,0]| = |{C G B^'° : Birth((:;) G [-t,0]}| < Y^^ti^) (6-3) 



We prove that this is finite with probability one by showing it has a finite mean. Indeed, 
reasoning as in the previous section. 



\ e / e 



(6.4) 



if a < 00. □ 
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7 Proof of Theorem |2.10| 



We prove that the hypothesis of Theorem 2.10| imply those of Theorem 4.6. The different 
parts of Theorem 2.10| follow by combining Theorem 4.6 and the space- width and time-length 



estimations of Theorem 5.27 



Existence and uniqueness In Theorem |5.27] (i) the condition a(/3) < 1 was shown to 



imply lack of backward oriented percolation. This, plus part (1) of Theorem 4^, proves (1) of 
Theorem ^TOj 



Exponential time convergence Inequality ( p.llD follows from ( |5.28D , ( |5.29D and ( |49| ) by 

choosing 6 = (2 — a)^^. □ 

Exponential space convergence In view of ( [4.21| ) it suffices to bound 

P(SW (A(Supp (/)))> d(Supp(/),A^)) < Yl IP(sW(A-'0) > d({x},A^)). (7.1) 

xeSupp(/) 

By ( ^.31| ) , this is bounded by 



E 



a;eSupp(/) 



ao(^) 
1 - a0) 



-(/3-/3)d({x},A'=) 



□ 



(7.2) 



Exponential mixing We shall use part (4) of Theorem We first show, in the next 
lemma, that a < oo implies the existence of an h such that there is non-oriented percolation in 
the interval (0, /i). 



Lemma 7.3 For all h > such that 

ah < 1, 

the probability that there is no (non- oriented) percolation in {0,h) is one. 



(7.4) 



Proof. Analogously to Section 5.1, we dominate the construction of the set ancestors of the 



non oriented percolation process by a multitype branching process. In this branching process, 
the mean number of ancestors of a contour 7 is 



fn{j,e) =hl{e '/"y}e 



-m 



(7.5) 
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As in Lemma 5.13, this branching process is subcritical ifa = a/i<l. □ 



We prove now From (|3.37D and (|3l^) : 



HK{fg) - ^a/ i^kq = mim)g{m)) - ^fivo) %(%)• (7.6) 

By ( 4.12| ) it is enough to bound 

p(a(Supp(/)) / A(Supp(5))) (7.7) 

where A(Supp (g)) has the same distribution as A(Supp (g)) but is independent of A(Supp (/)). 
This is bounded by 

p(a^'° 7^ AS^'°) < IP(SW(A^'°) + SW(A2^'°) > Ix-yj) (7.8) 

xeSupp{/), a;GSupp(/), 
i/GSupp (g) yeSupp (g) 



Using the following inequality which is valid for independent random variables Si and S: 



2, 



nsi + s2>i)<Y > j) >i-j) (7.9) 



and the exponential decay of ( ^.31| ) we get the decay stated in ( 2.13| ). □ 



Central Limit Theorem We apply the central limit theorem for stationary mixing random 
fields proven by Bolthausen (1982). Let = r^f- Let Aa be the sigma algebra generated by 
{Xj: : X e A}. Define 

ak,e{n) = sup{|P(^i n ^2) - P(^i)P(^2)| : G ^Ai, ^2 G ^A^, 

|Ai| < A:,|A2| <£, d(Ai,A2) > n}. (7.10) 

The simplified version of Bolthausen theorem stated in Remark 1 pag 1049 of his paper says 

that if there exists a 6 > such that ||Xa;||2+5 < 00 and 

00 

5^n'^-H«2,oo(n))^/(2+5) <^ (7.11) 

n=l 

then D < 00 and ( |2.14| ) holds. Hence, it suffices to show that 0:2,00 ('^) decays exponentially fast 
with n. We can write 

02,00 (n) = sup |/u(5i5'2) - W2I (7.12) 
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where the supremum is taken over the set of a € Z'^, gi in the set of indicator functions with 
support on Supp (/) U TaSupp (/) and (72 in the set of indicator functions with support in 

|j|TySupp(/) : y € Z'^ and \y - x\ > nVx G Supp(/)}. (7.13) 

By (|jD 



xgSupp (gi), 
j/GSupp (92) 

< 4(M2)2|Supp(/)| e""^''"' (7-14) 

|y|>n-2|Supp(/)| 

because |Supp(5'i)| < 2|Supp(/)| and || 

5i||oo — II52II00 — 1- Hence, O2,oo(^) decreases exponen- 
tially fast with n. □ 



8 Proof of Theorem |2.17| . Poisson approximation 



We define first a common probability space where all processes {r/^ :/?*</? < 00} can 
simultaneously be constructed. For each 7 G G = U^Gj let be a marked Poisson process on 
q£ -^^ The event points of this process are denoted {u,t,r,s), where u G M'^, t,r G M 
and s G M^. The coordinate t is interpreted as time while the coordinate r is later used to tune 
the rate of the projected process (u, t). The coordinate s — the mark — is an exponential random 
variable with mean one independent of everything (used later to determine the lifetime of the 
corresponding point/cylinder). Denote by P the product measure generated by (N^ : 7 G G), 
and by E the corresponding expectation. We identify the random counting measure with 
the corresponding discrete random subset of R"^^^ x M+. 

Fix a contour length j and an inverse temperature /3. By counting only those points in 
whose r coordinate is in [0, e"^*^!^!^-'''] we generate the (d + l)-dimensional marked process 

^e,i3 = marked Poisson process of rate e~^^^^^~^\ (8-1) 

The life of each point (n, t, s) G Ng^^ x R"^ is the interval [t, t + s]. 
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Define a family of marked point processes indexed by /3 and TxO, x G Z'^, 9 ^ G for Borel 
sets / C M by 

Nr^eAl) = N,,^(0(xe-^^Vrf,e-/3i/^/2) x /) (8.2) 
where for y = {yi, . . . , y^) G Z'^, O is the d-dimensional "rectangle" 

0{y, p) = [yi - p,yi + p] X . X [yd - p,yd + p] (8.3) 

Since the volume of 0{xe~^^^'^,e~^^^'^/2) is e~^^ , and Ng^^ has rate e~^^^^^~^\ the resulting 
process Nt^q^I) is a one dimensional marked Poisson process of rate e"^'^'. The marks are the 
independent exponentially distributed random variables of mean one, inherited from ^e,p- The 
point of this construction is that all these Poisson processes are constructed simultaneously as 
a function of the original d + 2-dimensional Poisson processes. 



Now we use the processes N^-^g^p to perform the graphical construction of Section 3.1. We 
call be the family of cylinders so obtained. Let be the free network of Section 3.2. and 
r/f the loss networks of (p.6|). As in (|3.8D and (3.37) these processes have invariant distributions 



pPp and pp respectively. 

Let y be a d-dimensional rectangle as in the statement of the theorem. Let 

Kp^^) = E ^oir.l) = E E l^B^^i^ = Life iC) 3 0} (8.4) 



as in ( |3.6| ). The super-label zero in the left hand side indicates we are dealing with the free 
process , while the sublabel zero in the right hand side indicates time zero. The family 
(M^^(y) : 7 S Gj) consists of \Gj \ independent Poisson random variables with mean: 

EM°^(y) = \V ■ e^l^l/'^l e-^l^l. (8.5) 

By ( |3.33|) ?/o(t) constructed with the cylinders in is pp distributed. Thus we can use 
77^(7) in the definition (|2l6l) of M^^p. By (p|) 



%^(7)<Co^(7). (8.6) 

Hence, M^^^V) < M^^V). 
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The joint construction also implies 



nMlp{V)-M,,p{V)>l) 



< 



< 



E IP(?o(^x7)-r?o^(r.7)>l) 
E k^o^(r.7)>l,%^(r.7)=0 

+ P(eo(Tx7) >2,r?^(r,.7) = 1 
From the construction, for any 9 £ G, 

{CoiO) >l,Vo{^) = 0^ C |c^:C^9CwithBasis(C) =eLife(C) 9 0,Af /0}. 
The probability of this last event is bounded by 

p(Co>)>i)ip(Eb?(^')>i)- 



Since ¥{^^{6) > 1) = 1 - exp(e-'3|^l) < e'^l^l, the rhs of Q is bounded above by 



On the other hand, 



P(?o^(e) > 2,r/o^(0) = l) < p(^o^(e) > 2) < ^ e- 



2/3|e| 



From (l8?7| ) -( Ml|) we get 

nM%{V) - M^,f3{V) > 1) < \V- e^l^l/'^l e-'^l^l (I7I «(/?) + ^ e'^l^l) 



~ e 



-2d/3 



(8.9) 



.10) 



^.11) 



^.12) 



To finish the proof of ( 2.18 ) we must show that ^ is close to a Poisson process. For 
I7I = j, let M^^oo count those points of the Poisson process N^^oo whose life contains the origin. 
The process M^^oo is a Poisson process in R"' of rate one. This is because the lifetimes are 
independent exponentials of mean one and for I7I = j, N^^oo is a Poisson process of rate one. 
The family {M^^oo : 7 £ Gj} inherits independence from {N^ : 7 € G}. For J C let 



J:a= {r G M'^ : ra G J+ [1/2,1/2]'^} C 



.13) 
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By definition, 



Then, as {V ■ e*/'^) : e^j/'^ C V, 

p(MO^(y) - M^^^iV) / O) < p('m^,oo (v^ \ [{V ■ e"^!'') : e^^' 



.14) 



> 



< 



Y\{{y ■ e^^''^) : 



.15) 



Inequality (|2l8|) follows from ( plD and (|8l5|) . 

Proposition 1.2 of Neveu (1977) and the comments below the statement of the Proposition 
say that the distribution on finite unions of d- dimensional finite- volume rectangles is enough to 
characterize a point process. Since the estimates ( p.l8| ) can be easily extended to finite unions 
of rectangles, the weak convergence follows. □ 
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